A note on late-time tails of spherical nonlinear waves 
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Abstract 

We consider the long-time behavior of small amplitude solutions of the semilinear wave equation 
0,| Ocp = (j)P in odd d > 5 spatial dimensions. We show that for the quadratic nonlinearity (p = 2) 

1-^ ■ the tail has an anomalously small amplitude and fast decay. The extension of the results to more 

general nonlinearities involving first derivatives is also discussed. 
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In a recent paper pj we studied the late-time tails of spherical waves propagating on 
even- dimensional Minkowski spacetime under the influence of a long range radial potential. 
Using perturbation methods we showed that in six and higher even dimensions there exist 
exceptional potentials which produce tails with anomalously small amplitudes and fast decay 
rates. The main purpose of this note is to show that anomalous tails exist also for nonlinear 
waves. 

We consider the semilinear wave equation with the power nonlinearity 

□0 = 0^ □ = 9,2-A, (1) 

in odd d > 5 spatial dimensions {p > 2 is an integer). We assume that initial data are small, 
smooth, spherically symmetric, and compactly supported 

0(O,r) = £/(r), dt(j){0,r) = eg{r). (2) 



It is well known that the corresponding solutions exist globally in time (see, for instance, [2]) 
so there arises a natural question: what is the asymptotic behavior of solutions for t ^ oo ? 
In the following we address this question using perturbation theory. Our starting point is 
the perturbation expansion 

= £0O + £Vl+^V2 + --, (3) 

where ecpo satisfies initial data ([2]) and all higher 0„ have zero initial data. Substituting this 
expansion into equation ([1]) we get the iterative scheme 

□00 = 0, (4a) 

□0P-1 = ro, (4b) 

□02p-2 = P0O Vp-1' ^^'^' (^^) 

which can be solved recursively. 

The general spherically symmetric solution of equation (4a) is given a superposition of 
outgoing and ingoing waves 

0o(t,r)=0-*(t,r) + 0g^^(t,r), (5) 

where 

' (2/-A:)! aW(«) _^„,,,_, 1 ^, , ,,+i (2/ - fc)! aW(t;) 



V'O IJ^'^J- Ml 2^ ,,., _ ,^,/ , _^xi_fc, </>0 l^'H- ,+i2^l i> 



^^^ k\{l - k)\ {v - uy-^ ' v-u V . / ^'^^ -fb^ k\{l - k)\ {v - uy-'' 

(6) 



and u = t — r,v = t + r are the retarded and advanced times, respectively (the superscript 
in round brackets denotes the k-th derivative). Here and in the following, it is convenient 
to use the positive integer index / defined hj d = 21 + 3 (recall that we consider only odd 
spatial dimensions d > 5). Note that for compactly supported initial data the generating 
function a{x) can be chosen to have compact support as well (this condition determines a{x) 
uniquely) . 

To solve equation (4b) we use the Duhamel formula for the solution of the inhomogeneous 
equation D^ = N{t,r) with zero initial data (see, e.g.. 



t t+r-T 



0(t,r) 



2ri+i 



dr / p''-'Pi{fi)N{T,p)dp, 



(7) 



\t-r-T\ 

where Pi{p) are Legendre polynomials of degree / and /i = {r'^ + p'^ — {t ~ r)^)/2rp (note 
that — 1 < /U < 1 within the integration range). Applying this formula to equation (4b) and 
using null coordinates r] = r — p and ^ = r + p we obtain 



<Pp-i{'t,r) 



1 



t+r t—r 



2i+3r 



i+i 



d^ {^-vy^'Piii^mv,odv: 



(8) 



where now p = {r^ + {^ — t)(t — ri))/r{S, — tj). If the initial data ([2]) vanish outside a ball 
of radius i?, then for t > r + i? we may drop the advanced part of 0o and interchange the 
order of integration in ^ to get 

oo t+r 



<t^P-i(,t,r) 



2^^ / ^^ / (^ - vY^'Piif^) w'iv, or d^ ^ 



(9) 



— oo t—r 

Substituting ([5]) into Qj and using the identity (see the appendix for the derivation) \ 

t+r 



di 



t-r 



i^-vr ^ ^2/ + i)!! [{t-vy-rr-' [ 1+3/2 

, ,,2(n-2)^r'+i / „ ,r? ,^ / 1 



t — 1] 



we get 



</'p-i(i,^) 



f{i+i)p-i 



Ii{p,0) + O[- 



(10) 



(11) 



^ As in [1| we use the notation x- :— 1, x- :— x ■ {x — 1) {x — {k — 1)), fc > 0. 



where 

2(/+i)(p-i)-i 

C{1,P) ■■= {-ly ^^^^^^,, [(/ + 1)(P - 1) - 2]^ (12) 



and 

+00 



Ii{p,q):= ha^\^)Y{a^^-^'\r^)Yd^. (13) 



The coefficient Ii{p,q) is the only trace of initial data. We point out that there is no loss 
of generality in putting the coefficient C{l,p) outside the square bracket in flTTl) because if 
C{l,p) = (which happens for p = 2), then the integrand over rj in iQ becomes a total 
derivative, hence the whole integral (Q vanishes for compactly supported initial data. 

It is not difficult to verify that generically 02p-2 and all higher-order iterates also decay as 
X/;(('+i)p-i^ thus 0p_i gives a good approximation of the full tail provided that e is sufficiently 
small. More precisely, for fixed r and t ^ oo we have 

0(t,r)?^£P0p_i(t,r), (14) 

up to an error of order 0(t~I('+^)P^^l)0(£:^P^^). We remark that for / = the series ([3]) 



Q 



was proven in [J] to converge for small enough e, thereby making the approximation ([T] 
rigorous. We have not been able to prove an analogous convergence result for / > 1, however 
from the practical point of view the asymptotic nature of the perturbation series is sufficient 
in using the approximation ( IT4|) to make quantitative physical predictions. 

The advantage of the approach presented above, in contrast to decay estimates in the 
form of inequalities, is that it makes easy to identify anomalous tails for which the amplitude 
of the leading order term in the perturbation expansion vanishes. In the case at hand, as 
mentioned above, this happens for the quadratic nonlinearity p = 2 since from (IT2l) the 
coefficient C{1, 2) vanishes for any / > 1. This implies that there is no tail in the first order 
or, put differently, the system of equations (4a) and (4b) satisfies Huygens' principle (note 
that this is not true in three spatial dimensions, i.e. for / = 0, [4]). 

In order to determine the tail for the quadratic nonlinearity we need to go to the second 
order. Applying the Duhamel formula (I7l) to equation (4c) we obtain 

t+r t—r 

Mt,r) = ^y:p^ I d^Ji^-vy^'Piifi)Mv.OMv,OdV- (15) 

\t-r\ -C 



To compute the asymptotic behavior of this expression near timehke infinity we need to 
find first the asymptotic expansion of (f)i{t,r) near null infinity {u = const and v -^ oo). 
Substituting (jH]) into (4b) we get 

h{u) „.. ,.,,.. . , . 22' ru 



(j)i{u, v) = free part + 



(V — u 



\2l- 



- + 0{l/v 



2l+2\ 



h{u) 



I 



/U 
[a^'\x)fdx. (16) 
-oo 



Plugging (IT6l) and ([6]) into flTSj) . interchanging the order of integration, and expanding in 
powers of 1/t, we get the leading order asymptotic behavior at timelike infinity 



02(t,r) = (-l)' 



yM 



1 



T,_i(l,2) + 



2/(2/ + 1) t3«+i 
Thus, for p = 2 the approximation flT^ should be replaced by 



0(t,r)^£302(t,r) + C(£^) 



(17) 



(18) 



where 02 is given by ( 1171) . 

Since quadratic nonlinearities are common in nonlinear perturbation analysis, the anoma- 
lous tail (TTSjl appears frequently in applications. We emphasize that these applications are 
not restricted to higher dimensions and might be physically relevant because some important 
equations in physics (we mean, in four dimensions) are equivalent to scalar wave equations 
in higher dimensions (actually, we first came across this phenomenon while studying the 
Yang-Mills equations in four dimensions [5|]). 

The analysis presented above can be readily generalized to incorporate nonlinearities 
involving derivatives, for example equations of the form 



00 = 0P(a9i0 + pdr^Y 



(19) 



where integers p and q satisfy p+q > 2, and a, [3 are constants. For this equation, proceeding 
along the same lines as in the derivation of the tail flTTj) . we obtain 



0p+5_i(t,r) = {a- (3Y 



C{l,p + q) 



Ii{p,q) + 



(20) 



There are several special cases when the tail (120|) vanishes and the decay is faster: 

• g = and p = 2. This case has been discussed above (see (iT711 ). 

• q = 1 and p > 1. In this case the coefficient Ii{p, 1) vanishes (since the integrand in 
(TT3|) is the total derivative and a(x) is compactly supported) and, instead of (!20l) . we 
have 



^p{t,r) 



^(i+i){p+i) 



Iiip+l,0) + O[- 



(21) 



where 
D{l,p) 



2(i+i)p-i 



(2/ + 



\^-a)'-^^M^X^-2^ 



p+l {l + l)p-l 



(22) 



Note that for p = 1 and /? = the coefficient D{l,p) vanishes. In this case there is 
no tail at the ffist order whatsoever, while at the second order, in analogy to fITTj) . we 
obtain 



l^2r,3l-2 



,(t,r) = (-l)'a^2 



3/ + 1 1 



2/(2/ + 1) t3« 



+2 



Ii{3,0) + O 



q = 2 and p = 0. In this case C{l,p + q) = and, instead of fl20l) . we have 



0i(t,r) = (-l)'a/3--— -(/ 



1 



XK2,0) + O(i 



(23) 



(24) 



(2/ + l)!!'^ ' ^' t2i+3 

thus for a/? = there is no ffist order tail, in analogy to the case q = and p = 2. At 
the second order we get 



02(t,r) 



-1 



d+i 



{a-PY 



)3/ 



Xz(0,3) + O 



(25) 



2/(2/ + 1) t3«+i 

We note that for / = 1 and y^ a y^ (3 ^ this case is exceptional in the sense that 
the ffist-order tail decays faster then the second-order tail. This is in fact a peculiar 
property of the nonlinearity of the form dt(f)dr4> in 5+1 dimensions, as {adt4> + l3dr4>) = 
a^ {dt(j)f + 2af3dt(j)dr(f) + (3^ (9,0) ^ 

Note also that from i^^ and fl25l) the nonlinearities {dt4>) and {dr4>) produce exactly 
the same tail, in agreement with the well-known fact that the equation n</) = {dt(f)) — 
{dr4>) is Huygensian. 



g > 1 and a = (3 ^ 0. In this case we have 

Eil,p,q) 



<Pp+q-l(t,r) 



l[l+l){p+q)+q-l 



Xi{p + q,Q) + 



(26) 



where 



2(«+l)(p+g-l)+g-l 
E(/,p,g) = (-l)^+V ^^^^^^,, {l + l)\{l + l){p + q-l) + q-2f. (27) 

The formula (126|) reduces to ( 12T1) if g = 1 and p = 1, and to f l2^ if g = 2 and p = 0. 

Finally, we wish to remark that all the above analytic predictions have been scrupulously 
verified numerically. 
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APPENDIX 

Here we derive the identity (fTOj) . Changing the integration variable from ^ to /i we get 

'"^^7^ = "^';'K, , Jdi^p^ii.) (i - ^^V" . (A.i) 

t-r -1 



Expanding the power and integrating using the identity 

,k_ ^ (2/ + l)/c 



d 



^•^^...S-.... ^^-^^-(¥)K^.^.l)i ^^^^' (A.2) 



we obtain 

t+r 






t-r 0^™ 



l + 2mj m\{2l + 2m + l)\\ \t - r] ^ 

(A.3) 
Expressing the sum in terms of the hypergeometric function we get (fTOj) . 
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